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The standard Black-Scholes theory of option pricing is extended to cope with underlying return 
fluctuations described by general probability distributions. A Langevin process and its related 
Fokker-Planck equation are devised to model the market stochastic dynamics, allowing us to write 
and formally solve the generalized Black-Scholes equation implied by dynamical hedging. A sys- 
tematic expansion around a non-perturbative starting point is then implemented, recovering the 
Matacz's conjectured option pricing expression. We perform an application of our formalism to the 
real stock market and find clear evidence that while past financial time series can be used to evaluate 
option prices before the expiry date with reasonable accuracy, the stochastic character of volatility 
is an essential ingredient that should necessarily be taken into account in analytical option price 
modeling. 
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There has been a great interest in the study of the 
stochastic dynamics of financial markets through ideas 
and techniques borrowed from the statistical physics con- 
text. A set of well-established phenomenological results, 
universally valid across global markets, yields the moti- 
vating ground for the search of relevant models 0. 0. [5} . 

A flurry of activity, in particular, has been related to 
the problem of option price valuation. Options are con- 
tracts which assure to its owner the right to negotiate 
(i.e, to sell or to buy) for an agreed value, an arbitrary 
financial asset (stocks of some company, for instance) at 
a future date. The writer of the option contract, on the 
other hand, is assumed to comply with the option owner's 
decision on the expiry date. Options are a crucial ele- 
ment in the modern markets, since they can be used, as 
convincingly shown by Black and Scholcs 4], to reduce 
portfolio risk. 

The Black-Scholes theory of option pricing is a closed 
analytical formulation where risk vanishes via the proce- 
dure of dynamical (also called "delta") hedging, applied 
to what one might call "log-normal efficient markets". 
Real markets, however, exhibit strong deviations of log- 
normality in the statistical fluctuations of stock index 
returns, and are only approximately efficient. 

Our aim in this letter is to introduce a theoretical 
framework for option pricing which is general enough 
to account for important features of real markets. We 
also perform empirical tests, taking the London market 
as the arena where theory and facts can be compared. 
More concretely, we report in this work observational re- 
sults concerning options of european style, based on the 
FTSE 100 index, denoted from now on as St H - 

To start with, we note, in fact, that the returns of the 
FTSE 100 index do not follow log- normal statistics for 
small time horizons. We have considered a financial time 
series of 242993 minutes (roughly, two years) ending on 
17th november, 2005. The probability distribution func- 
tion (pdf) p(uj) of the returns given by u> = ln(5* / St-i) , 



taken at one minute intervals is shown in Fig.l. We ver- 
ify that the Student t-distribution with three degrees of 
freedom conjectured by Borland |(J provides a good fit 
to the data. A slightly better fitting is produced if a 
smooth gaussian envelope is introduced to truncate the 
distribution at the far tails, viz: 
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co = Log[S(i)/S(i-l)] (x 10" 3 ) 

FIG. 1: The pdf of returns evaluated for one minute time 
horizons for the FTSE 100 index, along a period of about two 
years, ending on 11/17/2005. The solid line is the truncated 
Student t-distribution given by Eq. Q, with a — 2.3 x 10~ 4 
and b/a = 5. 



A time-dependent family of self-similar Tsallis distri- 
butions 0, p(u),t) oc (a 2 (t) +w 2 ) 1 /(i-?), with variance 
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a 2 (t) oc t 2 /^~ q \ was previously considered in option 
price modeling [6j. The time parameter t refers here to 
the time horizon of return fluctuations, while q ~ 1.5 
is indicated from empirical pdfs. Even though sugges- 
tive results were found by Borland for a set of option 
prices taken from the real market, the time dependence of 
the variance a 2 (t) disagrees with observations. It is well- 
established that the variance of return fluctuations grows 
linearly with the time horizon within a certain range of 
days, and that the return pdfs are not self-similar (there 
is a crossover to the log- normal profile). A further diffi- 
culty has to do with the matching between the theoreti- 
cal and observed volatilities. As reported in Ref. Q, for 
instance, the model volatility needed to reproduce the 
volatility smiles of options based on the S&P500 futures 
in june 2001 is set to 32.95%. However, as the analy- 
sis of the S&P500 futures time series reveals, the actual 
volatility for that period was considerably smaller, close 
to 14.5% 8]. 

The essential reason for the specific choice of a(t) in 
the above pdf is that it is related to a formerly known 
Langevin equation, which leads, on its turn, to a partial 
differential equation for option prices. An alternative 
and more general point of view - to be pursued here - 
is not to advance, a priori, any hypothesis on the form 
of the return pdfs, while still dealing with a Langevin 
description of return fluctuations. 

Le us assume, therefore, that at time t — the under- 
lying index is So and that its subsequent values are given 
by 

S t = S exp\jd + x(t)] , (2) 

where x{t) is a dynamic random variable described, at 
time t, by an arbitrary pdf p(x,t). We take x = at 
t = 0. Observe that it is a simple exercise to write, from 
Eq. J2J , the formal expression for the pdf of returns with 
time horizon t, in a statistically stationary regime. 

As an inverse problem, we are now interested to find a 
function / = f(x,t), so that p(x,t) be derived from the 
Langevin equation 

^ = f(x,t) V (t), (3) 

where 77(f) is a random gaussian field, defined by (r](t)} = 
and (r}(t)ri(t')) = 5{t - t'). Actually, it is not difficult 
to compute f(x,t) as a functional of p{x,t). We just 
write down the Fokker-Planck equation that is satisfied 
by p{x,t) 0, 

to obtain, from the direct integration of (J2J, 



The function f(x,t) is an important element in the op- 
tion pricing problem, within the dynamical hedging strat- 
egy of investment. Considering the writer of an op- 
tion contract who owns a number A of the underlying 
asset, dynamical hedging, a concept firstly introduced 
by Black and Scholes, consists of defining option prices 

V = V(S 7 t) as a function of the asset's value S and 
time t, so that A = dV/dS, and the "minimal" portfolio 
H(S, t) = SA(5, t) — V(S, t) is imposed to evolve accord- 
ing to the market risk-free interest rate r, as if H(S, t) 
were converted to cash and safely deposited into an in- 
terest earning bank account. Under these conditions, a 
straightfoward application of Ito's lemma leads to the 
generalized Black-Scholes equation, 

— V-rV + rS— V+ -f 2 S 2 — V = . (6) 

dt ds r ds 2 1 ' 

Exact solutions of this equation may be given in terms 
of statistical averages. For the case of call options which 
expire at time t*, we get 

V = exp[-r(i* - t)]<eOSexp(C) - £)(Sexp(C) - E)) , 

(7) 

where C — L dt(x + r — \f 2 ) is a random variable de- 
scribed by some pdf cf) = (/)((, T), with T = t* — t being 
the option's time to expire. One may check, by direct 
substitution, that J7J) solves, in fact, Eq. (JSJ). The put- 
call parity relation is also satisfied by our solution, as any 
bona fide option pricing formula should do. It is worth 
noting that the case where f(x,t) is constant corresponds 
to the standard Black-Scholes problem. 

In order to develop analytical expressions for V(S,t), 
let Z{\) = (exp(zAC)) be the characteristic function asso- 
ciated to the distribution (j>((, T). We may write, without 
loss of generality, Z(X) = Z a (X, 0)Z a (0, X)Zb(X), where 

Z a {Xi,X 2 ) = (expfiAx (jf dt(x +r) - A 2 ± jf dtf 2 )]). It 
is clear that Z a and Z\> can be computed, in principle, 
with the help of the cumulant expansion. However, a 
non-perturbative analysis may be readily addressed when 
one realizes that Z a (X, 0) is nothing but the characteris- 
tic funtion of p(x,T). The remaining problem, then, is to 
evaluate the corrections due to Z a (Q, X) and Zb(X) (they 
are in fact small in realistic cases). Up to first order in 
A, we have Z b (X) = 1 and Z a (0, A) = 1 - i| ff dt(f 2 ). 
Using ©, we get Z a (0,X) = 1 - i±([x 2 (t*) - x 2 {t)\) ~ 
exp{— i^([x 2 (t*)— x 2 (t)])}. Taking now, as a phenomeno- 
logical input, that the volatility depends linearly on the 
time horizon T, we write Z a (0, A) ~ exp[— z-^OqT], which 
implies that 

<XC,T) = p(C-rT + a 2 T/2,T) . (8) 

It is not difficult to verify that the above expression is 
exact when f(x,t) = ao, and leads to a n op tion pricing 
formula previously proposed by Matacz [lOj . which was 
based on heuristic arguments. However, the usefulness 
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we write, for a time series of temporal length T — Ne, 

1 f N ^ 

- / dtf 2 ~ -(xjy - so) - [N - >^ exp(>i + i - x t )} , 

(10) 

where Xi = x(ie + t). Substituting IjlUfl into J7J), we find, 
then, a pragmatical formula to work with numerical sam- 
ples. 



TABLE I: Listing of call option premiums taken on 02dec05 
( S = 5528.1, g = 0.81, a* = 6.1%), 06dec05 (S = 5538.8, 
g = 0.91, a* = 6.9%), and 09dec05 (5" = 5517.4, g = 0.94, 
a* = 7.1%). The risk-free interest rate is r = 4.5%. Options 
expired on 16dec05. The mean volatility measured between 
02dec05 and 09dec05 is a = 8.0%. 
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TABLE II: Listing of call option premiums taken on 19dec05 
(S = 5539.8, g = 0.78, a* = 5.9%), 03jan06 (5" = 5681.5, 
g = 0.83, a* = 6.3%), and 12jan06 (5 = 5735.1, g = 0.78, 
<7* = 5.9%). The risk-free interest rate is r = 4.5%. Options 
expired on 20jan06. The mean volatility measured between 
19dec05 and 12jan06 is a = 6.1%. 



of the approximation given by Eq. ijSJ is restricted to 
the cases where the far tails of p(x, t) decay faster than 
exp(— \x\), because of the exponential factor in Eq. 

0. 

In a more empirically oriented approach, we may at- 
tempt to compute V(S,i), as defined by Eq. JJJ, from 
real or Monte Carlo simulated financial time series, with- 
out having any information on the pdf p(x,t). Since we 
would not be entitled to use Eq. JSJ anymore, it is nec- 
essary to rewrite the expression J. dtf 2 that appears in 
the definition of £ in terms of known quantities. Observ- 
ing that Ito's lemma yields, from Eq. ©, 
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FIG. 2: Volatility smiles related to the option values taken 
on (a) 02dec05 and (b) 09dec05 (see Table I). The black and 
white dots refer to EP and MKT premiums, respectively. 



A stochastic process {xi] was generated from a high- 
frequency time series of the logarithms of FTSE 100 in- 
dex, m(S'j), consisting of 242990 minutes ending on 17th 
november, 2005 (it is essentially the same data used to 
establish the histogram depicted in Fig.l). The series was 
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TABLE III: Comparison between Black- Scholes (BS) and 
Monte-Carlo simulations based on the Student t-distribution 
(ST). We use the modeling parameters (S,a*,r) reported in 
Table I. Expiry date is 16dec05 as well. 

partitioned into 94 pieces of 2585 minutes each (roughly 
one financial week). In each one of these pieces, the 
market drift parameter fi was determined by the least 
squares method. Then, according to Eq. we define 
Sxi = x-i + i — Xi = hi(Si+i/ Si) — /i. We have found that 
improved results are obtained if extremely intense fluctu- 
ations are removed from the series, since they are likely 
to correspond to the market's reaction to unexpected 
events (we have checked this assumption in a number 
of cases) . Once the standard deviation for the whole se- 
ries is do = 2.6 x 1CP 4 , we removed from the series the 
fluctuations given by \8xi\ > 3 x 10~ 3 ~ 10cr (only 144 
elements were taken out from the Sxi series) . A sequence 
of 3000 samples, separated in time by one hour transla- 
tions, was used in each option price evaluation. It turns 
out that the samples are related to a period of reasonably 
well-behaved volatility (it fluctuates around some stable 
value). We also introduce a phenomenological factor g, 
so that trials with series of different mean volatilities are 
obtained in a simple way from the mapping Sxi — > gSxi. 

We have computed in this way option premiums based 
on the London FTSE 100 index (with e = 1 minute). The 
results are shown in Tables I and II, where market (MKT) 
and evaluated prices (EP) are compared. In the tables, 
NA stands for "not available" data, while X stands for 
prices we did not evaluate; the values between brackets 
correspond to instances where the time series extension 
was not large enough to get accurate predictions. The 
daily closing values of the FTSE 100 index, S, the g- 
factors, and the model and observed mean volatilities (ct* 
and a, respectively) are reported in the table captions. 
The numerical errors in the EP columns vary typically 
from 0.5% for the larger premiums to 10% for the smaller 
ones. When comparing the market and numerical results, 
one should keep in mind the existence of the usual bid-ask 
spread of option prices. 

As clearly shown in Fig. 2, we have been able to model 



volatility smiles which are in good agreement with the 
market ones. The remarkable feature of these results is 
that implied volatilities can be obtained from a single 
model volatility, a*, which is in fact close to the actual 
observed market volatility a. 

We have, along similar lines, evaluated option prices 
through Monte Carlo simulations. The stochastic process 
{Sxi} is now generated as independent events from the 
Student t-distribution with three-degrees of freedom. As 
shown in Table III, the Black-Scholes (BS) and Student t- 
distribution-based evaluations (ST) yield relatively close 
(and not good) answers for option prices. It is impor- 
tant to note the numerical results do not differ too much 
if Matacz's analytical option price formula is used with 
convolutions of the truncated Student t-distribution, Eq. 

, once the log- normal return pdf cores get large enough 
for time horizons of a few days. 

To summarize, we have provided strong evidence, from 
the empirical analysis of non-gaussian market data, that 
the combination of delta-hedging strategy and suitable 
Langevin modeling allows one to compute option premi- 
ums with reasonable confidence, from the use of Eqs. (0 
and On the other hand, we have found that both 

the Monte Carlo simulations and the delta-hedging an- 
alytical framework based on fat-tailed distributions and 
time-independent volatilities which are close to the ob- 
served averaged values, fail to predict real market option 
prices. Our results point out that efficient option pricing 
analytical tools have necessarily to deal with the stochas- 
tic nature of volatility fluctuations, a main distinctive 
aspect of financial time series. 

The author thanks Marco Moriconi for a critical read- 
ing of the manuscript. This work has been partially sup- 
ported by FAPERJ. 
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